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PROGRAM BRIEF 



^" Purpose - A program to be incorporated into FORTRAN II-D 

under Monitor I disk system as a relocatable library subroutine. 
Pseudo-random numbers from either a rectangular or a standard 
normal distribution can be generated. 



B. Method - The technique used is the power residue method. Each 
random number serves as the multiplicand for the next number. 
The formula is: 

10 

R. N. = K x R. N. ( mod. 10 ) 

n+1 n 



where K is a 10 digit constant multiplier defined in the sub- 
routine. Only the 10 low order digits of the products are re- 
tained as the random number. (See technique discussion in 
Statistical Forecasting for Inventory Control by Robert G. 
Brown, McGraw-Hill Book Co. , Inc. , New York, 1959, p. 164. ) 

C. Restrictions - This subroutine should only be used in programs 
in which the mantissa length of floating point numbers is 8 
(i. e. , F-8) . There is no restriction on the size of fixed point 
numbers. If the F size restriction can not be met, a separate 
FORTRAN II program can be written to punch a deck of random 
numbers which can be read as input into the main program as 
they are needed. 



D. Accuracy - The power residue method multiplies a starting 
value by a constant prime number to produce a product. It is 
only the low order half of this product which is used as the 
random number, i.e., P^ to P£q below. 

'» 

P l P 2 P 3 P 4 P 5 P 6 P 7 P 8 P 9 P 1Q P 11 P 12 P 13 P 14 P 15 P 16 P 17 P 18 P 19 P 

* random number 

The rightmost few digits of the product cycle in an easily 
detectable pattern and, hence, can not very well be called 
random. The "best" random digit is at Pjj, and it is this point 
from which the random number is taken. A cycle of 50 million 
10 digit numbers is expected before the sequence repeats. A 
more thorough discussion of the manipulation of this number 
is presented in the next section. 

- / — 
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E. Machine Configuration - Specifications for FORTRAN U-D 
must be met. The only special feature used in the subroutine 
is indirect addressing, which is a FORTRAN II requirement. 

F. Program Requirements - This subroutine uses 574 positions 
of core. It contains its own normalizing routine to store a 
rectangular (. 00000000 - . 99999999) or normal (generally 
-3. 5 to +3. 5) number in proper floating point form in FAC 

(a pseudo accumulator in which FORTRAN deposits numerical 
results, ) 



G. Source Language - This subroutine was written in SPS and 
assembled using SPS II- D. The condensed object deck is 
provided along with the source deck. 

H. Program Execution Time - Times in the subroutine to generate 
random numbers on a model I 1620 are as follows: rectangular 
distribution, 25 msec; normal distribution, 250 msec. During 
this time the argument is obtained, a number is generated and 
normalized, and the result is placed in FAC. Times will vary 
slightly depending upon the number of leading zeroes in the un- 
normalized result. 

I. Check-Out Status - This routine was tested and is being used by 
one installation at the time of submission of this program. 
Results of tests agree exactly with results from 1620-7. 0. 021, 

a random number generator subroutine written for FORTRAN 
W/ FORMAT. 



J. Sample Problem Running Time - The execution time of the 
FORTRAN object program is approximately 6 minutes on the 
1620 Model 1. 



K. Comments - This program and its documentation were written 
by an IBM employee. It was developed for a specific purpose 
and submitted for general distribution to interested parties in 
the hope that it might prove helpful to other members of the data 
processing community. The program and its documentation are 
essentially in the author's original form. IBM serves only as the 
distribution agency in supplying this program. Questions con- 
cerning the use of the program should be directed to the author's 
attention. 



DETAILED PROGRAM DESCRIPTION 



A. Rectangular Numbers - In FORTRAN a relocatable library sub- 
routine is called by mentioning its name and the argument upon 
which it is to operate, e. g. , SQRT(Y) or LOGF(2. 621). In 
1620 FORTRAN II-D the address of the argument is placed in 
the 5 positions of memory immediately preceding the location 
of the first instruction of the subroutine. The subroutine itself 
is placed, along with other subroutines called by the mainline 
program, at the end of the object code which FORTRAN gener- 
ates. (See pp. 133-134 of Monitor I Reference Manual, 
C26-5739-2, for rules on subroutine operations.) 

The power residue method of generating pseudo-random 
numbers requires a "seed", or starting value, to begin pro- 
ducing numbers. In this routine the "seed" is a 3 digit number, 
say, 231, in which case the user might write X=RAND(. 231). 
The argument is stored in floating point form as 2310000000. 
This is converted by the subroutine into 2310000001, which 
serves as a 10 digit multiplicand to be multiplied by 1977326743, 
a 10 digit prime number defined in the subroutine. The 10 
low order digits of the 20 digit product are retained. In this 
example they are 8307326743. 

This number is both the random number and the multipli- 
cand for the next desired number. It must now be placed in 
FAC in correct floating point form. The 2 low order digits 
are replaced by the exponent and the number is placed in FAC 
as 8307 326700, which is . 83073267. Should the low order 
half of the product read 0071622343, the subroutine would alter 
its form to 0071622300 and then begin to delete leading zeroes 
and adjust the exponent in the normalizing routine. Finally 
the normalized result 7162230002 is obtained in FAC. This 
is, of course, . 716223 x 10" 2 or . 00716223. 



B. Normal Numbers - Normally distributed numbers are generated 

from rectangular numbers in an interesting fashion. The Central 
Limit Theorem of statistics states that for a series of indepen- 
dent, identically distributed variates X}, X£, . . . x n with 
mean/£ and variance O r2 , the statistic 

S = X - JUL (where x" = £] x i ) 

«"//*" ~ 

or S = ^x - iijl 

ys- <r 

becomes normally distributed with mean and variance 1 as 
n — te-oo . For practical purposes very good results can be 
obtained with a modest sample size. 
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For a rectangularly distributed variable with limits of and 
1, the mean and variance of this variable can be shown to be 
. 5 and 1/12 respectively. If a sample of 12 such independent 
numbers is taken, this sample will have an expected value of 
4 and a variance of 1 since the variance of a sum of independent 
variables is equal to the sum of their variances. Since a 
rectangularly distributed random number has limits of and 1 
and since these numbers are theoretically independent, the 
statistic 

S = ^R. N. i - 12(.5) 

= Er. N. i - 6. 

can be expected to be approximately normally distributed with 
mean and variance 1. Therefor, it is seen that normal num- 
bers can be produced by summing 12 rectangular numbers and 
subtracting 6. 0, It is possible to produce both positive and 
negative numbers ranging from -6 to +6 with this technique. 
The theoretical normal distribution has limits of - oo and + oo , 
but, for all practical purposes, numbers outside the range of 
-3. 5 to +3. 5 are extremely rare. Tests of normality of a 
distribution of pseudo-random normal numbers have shown 
this generation technique to give quite satisfactory results. 



in. SUBROUTINE CALLING SEQUENCE 

The object deck for the subroutine is provided with name 
cards which allow the subroutine to be called RAND or RANDF 
in keeping with customary FORTRAN subroutine naming rules. 
Every time a statement mentioning RAND or RANDF is en- 
countered, a new number will be generated. A positive argu- 
ment signifies a rectangularly distributed aunber, e. g. , 
RAND(.179); a negative argument, a normally distributed 
number, e.g., RANDF(-. 179). 

As long as the argument in the calling statement does not 
change, the subroutine uses the last random number generated 
as the multiplicand for the next. This is done by saving the 3 digit 
argument in the subroutine and comparing it with the argument 
found each time the subroutine is called. Should the argument 
change, a new sequence is begun as defined in the preceding 
section of this write-up. (NOTE: . 179 and -. 179 are not the 
same arguments. ) 
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The example below will generate 50 rectangular numbers 
from one sequence, and then initiate a new sequence of 100 
normal numbers. 

DO 10 I = 1, 50 
10 A(I) = RANDF (.337) 

DO 20 I = 1, 100 
20 B(I) = RAND (-.113) 

The following example is frequently seen in simulation 
problems. Suppose copper anodes are normally distributed 
with a mean of 250 lbs. and a standard deviation of 3 lbs. 
Simulated sampling of such a distribution can be obtained by 
the FORTRAN expression 



WEIGHT = RAND (-. 707) * 3. + 250. 



UT. Rauo Subroutine: Program Flowchart 
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OPERATING INSTRUCTIONS 



A. Console Switch Settings - Same as FORTRAN II-D requirements. 

B. Loading Subroutine onto Disk - A condensed object deck with 
necessary control cards is provided for the user (deck 2). 
The structure of this deck is seen on p. 135 of the Monitor I 
Reference Manual. A monitor call and JOB card are also 
provided. The *DLOAD card (4th card) lists the subroutine 
name RAND in cc. 7-10. This can be changed to any name 
the user desires. The DIM number for this subroutine is 0026 
in cc. 17-20 of the same card. If the user has added or deleted 
library subroutines previously, this number may have to be 
changed accordingly. (See rules on p. 135 of Monitor I Reference 
Manual. ) The last card of the object deck contains the alternate 
name by which this subroutine may be called, RANDF, subject 

to the user's preference. The DIM number is also 26 in cc. 14-15. 
This may have to be changed as mentioned above. A complete 
listing of the object deck appears in another section of this 
write-up. To load the object deck merely depress the LOAD 
key at the card reader. 



. SAMPLE PROBLEM 

This sample problem generates and prints 10 rectangularly 
distributed random numbers and 10 normally distributed random 
numbers. It then generates 500 rectangular numbers and tal- 
lies a distribution of these numbers with . - . 09999999 as 

case 1, . 1 - . 19999999 as case 2, 9 - . 99999999 as 

case 10. Finally 500 normal numbers are generated and tal- 
lied with those below -3. as case 1, -3. to -2. 5000001 as 
case 2, ..... 2. 5 to 2. 9999999 as case 13, 3. or greater 
as case 14. At the conclusion of the sample problem "STOP" 
is printed at the console typewriter. The results of the sample 
problem appear in the Listings section which follows. 

The user should compile and execute the FORTRAN test 
problem on his own 1620-1311 system. This program tests the 
correct placement of the RAND subroutine onto the 1311 disk. 
(NOTE: the random number generator subroutine must be 
stored on disk before the program will execute successfully. ) 



o 



VII. PROGRAM LISTINGS 

The following program listings are presented: 

1) SPS source for RAND subroutine with control cards 

2) SPS condensed object deck with control cards 

3) FORTRAN program using RAND subroutine as sample 
problem (no data required) 

4) Results of sample problem 
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3400032007C 13600032007024902402 51 196 361 1300102 
ZZ JOB 5 
ZZSPS 
*L IBR 

•NAME RAND 

•ASSEMBLE RELOCATABLE 
•STORE RELOADABLE 
•10 NUMBER 0026 
• OUTPUT CARD 



1040 DSA RAND 

1070RAND SM ARGAD0 . 2 » .DECREMENT ADDRESS TO GET AT MANTISSA OF ARGUMENT 

1080 BNF C0MP. ARGADD. 11. BRANCH TO COMP IF MANTISSA IS POSITIVE 

1090 SM ARGADD. 5 ..DECREMENT ADDRESS TO GET AT 3 DIGIT ARGUMENT 

1110 SF ARGADD . . 6 .FLAG LOW ORDER DIGIT OF ARGUMENT SINCE NEGATIVE 

112 B C0MP612 

1130 00RG»-3 

114wC0MP SM ARGADD. 5. .DECREMENT ADDRESS TO GET AT 3 DIGIT ARGUMENT 

1150 C STORE .ARGADD. 11. COMPARE OLD ARGUMENT WITH NEW 

1160 BE SUMSET ... IF EQUAL CONTINUE CHAIN OF R.N./S 

1181, TF STORE .ARGADD. 11 .SINCE UNEQUAL PUT NEW ARGUMENT INTO STORE 

1190 T- R NUM. CONST . .LOAD ZEROES AND 1 INTO RNUM TO FORM MULTIPLICAND 

1200 T- RNUM-7. STORE. .LOAD 3 DIGIT ARGUMENT UNTO RNUM 

203C CF RNUM-7. ♦ .REMOVE FLAG IF IT EXISTS 

2040SUMSETTF SUM, NE GS I X . .SET SUMMING AREA TO -6 

2055 TFM CTR. 12.10, SET COUNTER TO 12 . 

2 J7JCALC M RNUM. KL . .DEVELOP PRODUCT 

2C80 SF 0009C,,,FLAG HIGH ORDER END OF LOW ORDER HALF OF PRODUCT 

2093 TF RNUM. CO 99. .MOVE RANDOM NUMBER TO RNUM 

2110 BNF SINGLE .STORE. .3RANCH IF RECTANGULAR DIST. NUMBER DESIRED 

2120 A SUM, 000 )9 . » ADD NJMBER TO SUM 

2130 SV C TR » 1 . 1 Z , DECREMENT COUNTER BY 1 

2140 BNZ CALC. ..RETURN IF 12 NUMBERS NOT YET SUMMED 

2150 TF FAC-2 .SJM-2 . .MOVE OUT 9 HIGH ORDER DIGITS 

2160 BNF S I NGLE & 24 .SUM . .BRANCH IF RANDOM NORMAL NUMBER IS POSITIVE 

2170 SF FAC-2 ,. .OTHERWISE SET FLAG FOR NEGATIVITY 

2180 6 SIN5LE624 

2190 0ORG»-3 

3030SINGLET*-" FAC-2 .00097 . .LOAD RECT. DIST, NUMBER INTO FAC 

3040 TOM FAC- 10. 0,11 

3050 TFM FAC. 0,10, SET EXPONENT EQUAL TO 00 

3C8CNORM TO SI GN, FAC-2,, THIS IS NORMALIZING ROUTINE USED BY BOTH 

3090 8D ADDFXP.FAC-10 

3100 CF FAC-2 

3110DIGTCKBD FlGSET,FAC-9 

3120 TO SAVE, FAC-2 

3130 SM FAC. 1,10 

3140 TD FAC-2, RKMK 

3150 TR FAC-9,FAC-8 

3160 TO FAC-3.SAVE 

3170 TDM FAC-2, 

3180 B DIGTCK 
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RAND S oaAGO Ti n£ (C on r. ) 
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DS 
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OS 
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OS 
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DS 
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1060ARGAODDS 


.RAND-ltFIELD TO CONTAIN ADDRESS OF ARGUMENT OF 


4220 


0END1 
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RAND 3 ubroutine — Object 
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FORTRAN E-D SAMPLE i-'RCKLm 



SAMPLE PROBLEM - TEST OF RAND 01 

DIMENSION KASE«14n 2 

PR IflT 100 03 

100 FORMAT *45HSOME RECTANGULARLY DISTRIBUTED RANDOM NUMBERS/n 04 
DO 1 MAX # 1,10 05 
X # RAND«.431n 06 

1 PRINT 101, X 07 

101 FORMAT *F11.8n 08 
PRINT 102 09 

102 FORMAT %/MOHSOME NORMALLY DISTRIBUTED RANDOM NUMBERS/o 10 
DO 2 MAX # 1,10 11 
X # RAND*-.483n 12 

2 PRINT 101, X 13 
PRINT 103 U 

103 FORMAT %///39HDISTRIBUTION OF 500 RECTANGULAR NUM3ER5//n 15 
PRINT 104 1 6 

104 FORMAT %13HCASE TALLY/a 17 
DO 19 I # 1,10 18 

19 KASE36 In #0 19 
DO 20 MAX # 1,500 20 
X # RAND*.733n 21 
1 # X*10,0 & 1.0 22 

20 KASESIa # KASE%In & 1 . 23 
DO 21 I * 1,10 24 

21 PRINT 105. I, (CASE* In 25 

105 FORMAT »I3,I9o 26 
DO 29 I # 1,14 27 

29 KASE%In #0 28 
DO 30 MAX # 1,500 ?9 
X # RAND%-.653n 30 
IF »X&3.0n 31,32*32 31 

31 KASESla * KASEfcla & 1 32 
GO. TO 30 33 

32 IF %X-3.Qn 34,35*35 34 

35 KASE»14n * KASE*l4n & 1 35 
GO TO 30 36 

34 I ft 2.0*X & 8.0 37 

K.AS £% In * KASEVIa & 1 38 

30 CONTINUE 39 
PRINT 106 40 

106 FORMAT »///34HDISTRIBUT ION OF 500 NORMAL NUMBERS//D 41 
PRINT 104 42 
DO 36 I # 1,14 43 

36 PRINT 105* I, KASESin 44 
STOP 45 
END 46 
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Assembly 



L ISTIN© 
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dsa 


01070 


RAND 


SM 


01080 




RNF 
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SF 
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COM? 


C . 
DORG 
SM 


01150 




C 


0116 




BE 


01180 




TF 


01190 




TF 


01200 




TF 


02030 
0?0i»0 
02 050 
02070 
07080 


surtSET 

CALC 


CF' 

TF 

TFM 

M 

SF 


07090 




TF 


02110 




RNF 


07120 
02130 
0714 




A 
SM 

mz 


07150 




TF 


07160 




8MF 


02170 




SF 


07180 
07130 
03030 


B 

00 RG 
SINGLE* TF 


03C4 
3 05 
050«0 


iior.;i 


TWt 
TFM 
TO 


G309G 
03100 




RD 
CF 



RAND 00001* 00005 3000G 

ARGA DO , 2 , , DEC REMENT ADDRFSS TO GET AT MANTISSA OF ARGUMENT 

00006 12 00005 7J0002 
COMP,ARGADD, 11, BRANCH TO COMP IF MANTISSA IS POSITIVE 

0001C i»it 00062 OOODf 
ARGADD, 5,, DECREMENT ADDRESS TO GET AT 3 Din IT AHRWEHT 

00030 12 00005 B0005 
ARGA00,,6,FLA« LOW ORDER DIGIT OF ARGUMENT SINCE NEOAT IVE 

0001*2 32 00005 00000 

C0MP*I2 0005% 49 00071* 00000 

*-3 00062 

A RGADD, 5 , ,PEC REMENT ADDRESS TO GET AT 3 DIGIT ARGUMENT 

1 00062 12 00005 S00G5 

STORE *ARGADO, 11, COMPARE OLD ARGUMENT WITH NEW 

0007U 55 00367 00005 
SUMSET,,,IF EQUAL CONTINUE CHAIN OF R.N./S 

00086 1*6 0011*6 01200 
STORE ,AAGADD, 11 .SINCE UNEQUAL PUT NEW ARGUMENT INTO STORE 

00098 56" 0036 7 00005" 
RNUM, CONST,, LOAD 2EROE.S AND 1 INTO RNUM TO FORM MULTIPLICAND 

00110 55 00533 0051*1 
RNUM-7,ST0RE,,L0AD 3 DIGIT ARGUMENT UNTO RNUM 

00122 55 00526 00367 

RNUM- 7 , , , REMOVE FLAG IF IT EXISTS 00131* 33 00526 00000 

SUM,NEQS IX , ,SET SUMMING AREA TO -6 00146 25 00573 0056 2 

CTR,1?#10,SET COUNTER TO 17 00158 16 003G9 00012 

RNUM;K,, DEVELOP PRODUCT 00170 53 00533 00551 

00090,,, FLAG HIGH ORDER END OF LOW ORDER HALF OF PRODUCT 

00182 32 00090 00000 

RNUM, 00099,, MOVE RANDOM NUMBER TO RNUM 

00194 ?6 00533 0009? 
SINGLE, STORE,, BRANCH 'IF RECTANGULAR DIST. NUMBER DESIRE" 

0070S £5 00798 0035 7 
0021S 5l 00573 0002" 
00230 1? 0036 9 OOOfll 

00242 57 00170 01200 

00254 26 02490 00571 
S»NGLE*24,$UM,, BRANCH IF RANDOM NORMAL NUMBER IS POSITIVE 

00266 £5 00522 00573 

FAC- 2,,, OTHERWISE SET FLAG FOR NEGATIVITY 

00278 32 02490 00000 

SINGLE+24 002T0 49*'00322 00000 

*-3 00293 
FAC-2, 00097,, LOAD RECT. OUT. NUMBER INTO FAC 

• 00298 26 02490 G0097 

FAC-10,0,11 00310 15 074S2 00003 

FAC, 0,10, SET EXPONENT EPUAL TO 00 00322 16 02492 00000 

SIGH, FAC-2,, THIS IS NORMALIZING ROUT ME USED DY T!OTI! 

00334 25 00493 07490 

APDEXP,FAC-10 00345 43 004C2 07402 

FAC-2 00353 33 07490 00000 



SUM, 00099,, ADD NUMBER TO SUM 
CTR,1, 10, DECREMENT COUNTER CY 1 
CALC,,, RETURN IF 12 NUMBERS MOT YET SUMMED 

FAC-2,SU«-2,*MOVE OUT 9 HIGH ORDER DIGITS 



..liUO 


r»i".*"C!'. 


"."> 


rLwsrT,r.*>c-o 


C0370 


•i 3 


GG4 CO 


fT>403 






in ' 


SAVE, FAC- 2. 


G03-? 


55 


Go4"i. 


i.7i|i0 






3:1 


r,\c,i,io 


003? 4 


12 


o^c 


oouSi 


v. ji-'i G 




f 


r*c-?,r.KMK 


CC4'.'i 




C'iCG 


or'*- 


j a j " 




rn 


FAC-9, r AC-C 


•004 lu 


31 


0'»4D3 


G'4 04 






Ti 


FAC- 3, SAVE 


GC430 


?5 


G74 0" 


0&4 C4 






r p;i 


FAC- 2,0 


0044 2 


15 


G'4TO 


00000 


ujioO 






niOTcr. 


00454 


59 


Gu37u 


GCCGO 






Dona 


i 


G04G2 








04 G3G 


ATTEaP 


AM 


FAC, -1,10 


0040 2 


11 


U'432 


OOCwl 


04 C4 C 




rr 


rAC-.7,FAC-3 


004 74 


20 


024?G 


G74S9 


u4 G50 


rL".o r " r 


cr 


FAC -9 


004 06 


32 


0?4o3 


OOGGb 


04 Go G 




'.;f 


*♦ 24, SIGN 


004 9 S 


45 


0G52 2 


004 33 


C.7C 




sr 


FAC-2 


G0510 


32 


07430 


GGUOO 


04 Cou 




rn 




C0522 


42 


GOCGO 


OOOGG 


o4 G'"0 




none 


*- n 


00524 












~>r. 


3,0,DIGTCK-3 


00307 


00003 uOO 


04 1?C 


:.:vr 


PS 


10 


00533 


00010 




G4130 


Cn.JST 

K 


"»C 


3,1 


00541 


00000 C.O00C001 


u4 14 


x 


10,1177326743 


00551 


00010 1977326743 


0415 


ix 


PC 


11,-60000000000 


00562 


0001B SOOOQOQQOO" 


04 luO 


sir 


ns 


11 


005 73 


00011 




G417 


ZTP 


ns 


2,nir,TCK-l 


00 36 9 


00007 




u4100 




DS 


1,FLGSFT*7 


00493 


00001 




04 KO 


SAVE 


OS 


1,FLGSFT*8 


004 94 


00001 




04 700 




DC 


1,0,FLGSET*9 


00495 


00001 * 




04 710 


FAC 


OS 


,249? 


07492 


00000 




0106 


ARGA^D 


DS 


,RANP-1, FIELD TO CONTAIN 


ADDRESS OF ARGUMFNT OF SUPR 








00005 


00000 




G4770 




TIJD 


1 


00001 
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Correction 
7.0.057 



Attached is a Comments Sheet 
\7hich should be included in 
the program writeup. 



/ 



v 



A user has recently pointed out that the following sequence of 
FORTRAN Statements will cause a check stop: 



Y = RAND ( X ) 
X = X + 1. 

The second statement will fail to execute properly because the 
generator subroutine sets a flag withing the mantissa of the argument, 
This occurs only when X is negative, i. e. , random normal numbers 
are being generated. 

It is suggested that the use of arithmetic operations be avoided . 
\ on the argument of the RAND subroutine. If such arithmetic is necessary 
(although the 'seed' is generally initialized and not altered during most 
runs) the following code may be used: 

X = Z 

Y = RAND ( X ) 
Z = Z + 1. 

X = Z 



